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We report on magneto-oscillations in differential resistivity of a two-dimensional electron system subject to 
intense microwave radiation. The period of these oscillations is determined not only by microwave frequency 
but also by its intensity. A theoretical model based on quantum kinetics at high microwave power captures all 
important characteristics of this phenomenon which is strongly nonlinear in microwave intensity. Our results 
demonstrate a crucial role of the multi-photon processes near the cyclotron resonance and its harmonics in the 
presence of strong dc electric field and offer a unique way to reliably determine the intensity of microwaves 
acting on electrons. 
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Over the past decade an array of remarkable effects 
was discovered in very high Landau levels of high mobil- 
ity two-dimensional electron systems (2DES). Among these 
are microwave-induced resistance oscillations (MIRO) JH, 
phonon-induced resistance oscillations J2], Hall field-induced 
resistance oscillations 

Hi, 

zero-resistance states [5], and 
zero-differential resistance states Theories of magneto- 
resistance oscillations are based on the quantum kinetic de- 
scription and consider two main mechanisms: 1) the "dis- 
placement" mechanism originating from modification of im- 
purity scattering by microwave (ac) or dc electric fields 
and 2) the "inelastic" mechanism stepping from the formation 
of the non-equilibrium energy distribution I1 1011 . MIRO are 
controlled by the ratio of microwave frequency u> = 2wf to 
the cyclotron frequency uj c = eB/m*, e ac = u>/oj c - Hall 
field induced oscillations, which appear in differential resis- 
tivity, are governed by edc = e£dc(2R c )/h^Uc, where £dc is 
the Hall field, 2i? c = 2vf/u> c is the cyclotron diameter, and 
vp is the Fermi velocity. Finally, in a 2DES subject to both 
ac and dc fields, the resulting oscillations were found to de- 
pend on simple combinations of ac and dc parameters, i.e. 
e ac ± edc lllUll2ll . Such dependence indicates that the rele- 
vant inter-Landau level scattering processes involved a single 
photon. 

Importance of processes involving multiple microwave 
quanta was suggested by numerous experiments 111311 report- 
ing MIRO-like features in the vicinity of fractional values 
of e ac . The most prominent series of, so called, fractional 
MIRO occurs near subharmonics of the cyclotron resonance, 
e ac = 1/2,1/3,1/4,.... Theoretical proposals considered 
both multi-photon Il4lll5ll and single-photon lfl5l[l6ll mecha- 
nisms to explain the response at fractional e ac . On the other 
hand, relevance of multi-photon processes near the cyclotron 
resonance and its harmonics was not considered either exper- 
imentally or theoretically. 

In this paper we report on another class of magnetoresis- 
tance oscillations in a high-mobility 2DES exposed to high- 



power microwave radiation and strong dc electric field. These 
oscillations are manifested by a series of multiple maxima and 
minima all occurring in the close proximity to the cyclotron 
resonance and its harmonics. Furthermore, the period and the 
phase of these oscillations depend not only on e ac and edc but 
also on the microwave intensity. This characteristic sensitivity 
to microwave intensity sets apart this strongly non-linear phe- 
nomenon from all previously reported resistance oscillations. 

To explain our experimental findings we propose a theoret- 
ical model based on quantum kinetics which captures all im- 
portant characteristics of the phenomenon. As we will show, 
this unusual effect owes to the quantum oscillations in the den- 
sity of states and a crucial role played by multi-photon pro- 
cesses. In the presence of radiation the electron states are split 
into Floquet subbands separated by fku in a similar way as 
quasi-particle states split in the course of photon-assisted tun- 
neling across a Josephson junction [laj. The scattering rate 
off disorder is then controlled by the overlap of such subbands 
leading to oscillatory behavior in differential magnetoresis- 
tance. 

Our experiment was performed on a 100 /im-wide Hall 
bar etched from a symmetrically doped GaAs/AlGaAs quan- 
tum well. After brief low-temperature illumination with vis- 
ible light, density and mobility were n e ~ 3.8 x 10 11 cm~ 2 
and \i ~ 1.3 x 10 7 cm 2 /Vs, respectively. All the data were 
recorded under continuous irradiation by / = 27 GHz mi- 
crowaves in a Faraday geometry at T ~ 1.5 K. Differential 
resistivity, r(J) = dV/dl, was measured using a quasi-dc (a 
few Hertz) lock-in technique. 

As shown in Fig.[T](a), the magnetoresistivity measured 
under microwave irradiation in the absence of dc field ex- 
hibits sharp MIRO features around the cyclotron resonance 
(e ac = 1) and its second harmonic (e ac = 2). In both cases we 
observe exactly one minimum and one maximum positioned 
roughly symmetrically about e ac = 1,2 (cf. vertical lines). 

In Fig.[T](b) we present differential magnetoresistivity r{B) 
measured at dc currents / from 32 pA to 78 /iA (in step of 4 
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FIG. 1: [color online] (a) Magnetoresistivity p(B). (b) Differential 
resistivity r(B) measured at dc currents from / = 32 to 78 fiA in 
step of 4 fiA. All data are acquired under irradiation by / = 27 
GHz microwaves at T = 1.5 K. Vertical lines mark positions of the 
cyclotron resonance and its second harmonic. 

/iA) at the same radiation frequency and intensity. Remark- 
ably, the data readily reveal multiple maxima and minima oc- 
curring in the proximity to the cyclotron resonance. For ex- 
ample, the data obtained at I = 64 /iA exhibit three max- 
ima (cf.,t) and three minima (cf.,j). With increasing / all 
oscillations move to higher magnetic fields in a regular fash- 
ion. While not so pronounced, similar behavior is observed 
near e ac = 2. These findings indicate that these magneto- 
oscillations are qualitatively different from all reported previ- 
ously. As we show below, the phenomenon owes to a non- 
trivial role played by multi-photon processes at integer e ac in 
the regime of strong microwave and dc electric fields. 

To interpret the experimental data, the dynamical screen- 
ing of microwaves by 2DES El 1711 has to be considered. 
Our theoretical results are expressed through the dimension- 
less microwave power V u which includes this effect explicitly, 



V ± = 



(2ei? c 



eocft 2 Tom + (UJ ± UJ c f 



(1) 



Here, S — £oc£ ac is the energy flux carried by the microwave 
radiation, r cm = eo(y / e + l)m*c/2ne 2 is the electromagnetic 
damping time, and e = 12.8 is the dielectric constant of GaAs. 
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FIG. 2: [color online] Offset of the resistivity peak near the cyclotron 
resonance, S — e^ c — 1 (open circles) vs. edc- Linear fit (solid line) 
performed for edc > 0.6 (solid circles) crosses 5 = at edc = — 1/8 
(dashed line) as prescribed by Eq. l(3j. 

In our 2DES, r cm ~ 3.8 ps and w c r cm ~ 1. For that reason 
the parameter has a weak magnetic field dependence near 
the classical resonance uj c = 2nf, which is the focus of the 
present study. At the same time, we can still retain only the 
active component V~ of the incoming radiation (below the 
"— " sign is omitted). At uj w ui c , V u — 1 corresponds to 
a microwave flux of one photon passing through the area of 
\ 2 H = Hc/eB during the time r om . 

We now present the main result of our calculations which 
is the expression for the differential resistivity Sr = r — r £>'■ 



Sr _ (4A) 2 r tr 
r D itTn 



COS 



2ire dc J U 



^vVu sin 7re a 



Vui sin 27redc cos 7re ac J\ 4 yl V u smite 



(2) 



Here, J n (x) is the Bessel function of n-th order, r tr and r n are 
transport scattering and backscattering times, ro is the Drude 
resistivity, and A = exp[— 7r/w c r q ] is the Dingle factor. In our 
2DES, the quantum scattering time r q ~ 20 ps, which corre- 
sponds to A 2 ~ 0.15 at the cyclotron resonance. Equation (0 
holds at high temperatures, T > ftw^/VZ and strong direct 



currents, e& 



> 



We now discuss the evolution of the oscillation maxima po- 
sitions with dc field. At Vui ^ 1, Eq. © reduces to 



Sr 



cx cos 



AyPu sin(7re ac ) - 7r/4 - 27re dc 



Thus the position of the maximum closest to the cyclotron 
resonance is given by e+ = 1 + 5, where S depends linearly 
on e dc : 



£(edc 



1 



2VK 



1 



(3) 



This expression allows us to extract the value of which 
cannot be reliably determined in conventional, e.g. MIRO, 
experiments. In Fig. [2] we present experimental evolution of 
the phase S = e+ — 1 with increasing edc and observe linear 
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FIG. 3: [color online] (a) Differential resistivity r vs. e ac measured 
at dc current I — 54 under irradiation by microwaves, (b) 
Correction to the differential resistivity calculated with Eq. 10 us- 
ing Vuj — 4.2 and I — 54 fiA. Vertical line marks position of the 
cyclotron resonance. 



dependence for edc > 0.6. The linear fit (solid line) gener- 
ates — 4.2 and passes through (<5, edc) = (0,-1/8), as 
prescribed by Eq. ©. 

Having extracted V u and obtained excellent agreement for 
the evolution of the peak position with ed c , we can now com- 
pare the experimental and theoretical curves directly. Fig.[3ja) 
shows our data obtained in microwave-irradiated 2DES carry- 
ing direct current I = 54 pA, which corresponds to ed c — 
1.65 at e ac = 1, In Fig.|3lb) we present the results of calcu- 
lations using Eq. (O with V u = 4.2 extracted earlier. We ob- 
serve that the theoretical curve recreates experimental data, in- 
cluding the number of measured oscillations, reasonably well. 
The theory also captures the amplitude of the observed os- 
cillations. Using Eq. (O with t^/tt, w 0.18 measured in 
recent non-linear dc transport experiments |4[ we arrive at 
Sr/ro — 0.1. This value agrees well with the experimen- 
tal data in Fig. a) and further reinforces the validity of our 
model as well as the procedure for extracting V u - 

The remaining part of the paper is devoted to the deriva- 
tion of the main theoretical result, Eq. (|2). Our derivation 
is based on the quantum kinetic equation of electrons in the 
presence of in-plane electric fields in the limit of overlap- 
ping Landau levels. We find the distribution function f e - v 
which depends on the direction of the momentum p = ppn^ 
of an electron at the Fermi surface in the quasi-classical ap- 
proximation, Ep ^> u> c . The current density is given by 
j = 2evF J cos(pv(e)f £;v ded(p/2ir. 

The distribution function satisfies the kinetic equation [8, 
II 



^cd v fe;ip — St{/} e;(/ 



(4) 



with the collision term given in the Floquet representation IU9I1 



N J 



-Nhw)-f v (e)}. (5) 



Equation (© accounts for scattering events whereby the par- 
ticle changes its momentum by Aq = p-p (n v — n v ' ) and 
absorbs (N > 0) or emits (N < 0) \N\ photons. The ki- 
netic energy change includes the work done by the dc elec- 
tric field as the result of shift of the guiding center, W vv > — 
eER c (sin ip — sin ip'), and the energy of absorbed or emitted 
photons. The scattering rate is 



(AT) 
T vv > 



(6) 



where 
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the density of states is v(e) = i/q [1 — 2A cos(27re/fiu; c )] with 
isq = m/h 2 ir being the DoS per spin. In the classical limit 
of constant DoS the normalization condition J2n = 
1/t w / holds, as the radiation changes electron energy dis- 
tribution while not affecting the total number of particles scat- 
tered into a given angular range. We look for the solution of 
Eq. © in the form / = fa + A/ cl + AA (1) /, where fa is the 
Fermi-Dirac distribution function, fa + A/ C i is the solution 
of the classical kinetic equation with A = 0, and 



A (1) /(£,^) 



dfaje) 2-ke 
A\ cos ip — cos - — 

as njjjr 



(8) 



We substitute Eq. © to Eqs. ©, © and © and neglect 
terms containing derivatives of the Fermi distribution of the 
second and higher orders since the contributions of such terms 
vanish upon the integration over the energy. We obtain 



K {N ) 



N 



(9) 



ip V > 



Tipcp' 



{W vip i + Nhu) cos 



2tt {W w , + Nhw) 



The angular integrations in Eq. (O can be carried out using 
the stationary phase approximation. At e^ c > maxjv^J, 1} 
the main contribution to the integral comes from narrow in- 
tervals centered at ip = ±7r/2, (p' — z F 7r /2- This corre- 
sponds to backscattering of electrons in strong dc electric 
field. To calculate A\, we substitute from Eq. (|7]l to 
Eq. (O and sum over N. The summation results in a set 
of (^-functions 6(wt ± 2irui/u> c + 2irk) and their derivatives 
(k = 0, ±1, ±2, . . . ). The remaining time integration repro- 
duces the final result of this calculation presented in Eq. (|2]). 

Without microwave radiation or in the weak power limit, 
Eq. (O reduces to the result of Refs. Jgt] or | L2], respectively. 
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At stronger microwave powers, the wave functions are of 
the Floquet type[ 19] and characterized by an additional time- 
dependent phase factor. According to Eq. (0, the number of 
photons maximizing a^S is N m — \fVZ, since in this case 
the phase factor oscillates in-phase with factor exp(iNujt) and 
the most efficient scattering off disorder is accompanied by 
absorption or emission of N m quanta of energy hut, Eq. ©. 

In summary, we reported on experimental and theoretical 
studies of magnetotransport properties of high-mobility two- 
dimensional electron systems driven away from equilibrium 
by intense microwave radiation and strong dc electric fields. 
We observed a new class of magnetoresistance oscillations 
near the cyclotron resonance and its harmonics, with positions 
of peaks and dips strongly dependent on microwave intensity. 
We proposed a theoretical model based on the quantum kinetic 
equation taking full account of multi-photon processes. Our 
calculations capture all important characteristics of the phe- 
nomenon — the period, the phase, and the amplitude of the 
oscillations are all in excellent agreement with experimental 
observations. Taken together, these results demonstrate the 
crucial role of the multi-photon processes near the cyclotron 
resonance and its harmonics in the presence of strong dc elec- 
tric field and offer a unique way to reliably determine the mi- 
crowave intensity seen by 2D electrons. 

We thank I. A. Dmitri ev and B. I. Shklovskii for useful dis- 
cussions and remarks. The work at Minnesota was supported 
by NSF Grant No. DMR-0548014. M.K. acknowledges sup- 
port by the US DOE under contract number DE-AC02-98 CH 
10886 and BNL grant 08-002. M.V. acknowledges the Donors 
of the American Chemical Society Petroleum Research Fund 
for partial support of this research. 



[1] M. A. Zudov, R. R. Du, J. A. Simmons, and J. L. Reno, Phys. 
Rev. B 64, 201311(R) (2001); P. D. Ye at al, Appl. Phys. Lett. 
79, 2193 (2001); R. G. Mani atal, Phys. Rev. Lett. 92, 146801 

(2004) ; J. H. Smet at al, ibid 95, 116804 (2005); A. T. Hatke, 
M. A. Zudov, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 
102, 066804 (2009). 

[2] M. A. Zudov at al, Phys. Rev. Lett. 86, 3614 (2001); A. A. 
Bykov, A. K. Kalagin, and A. K. Bakarov, JETP Lett. 81, 523 

(2005) ; A. T. Hatke, M. A. Zudov, L. N. Pfeiffer, and K. W. 
West, Phys. Rev. Lett. 102, 086808 (2009). 

[3] C. L. Yang, J. Zhang, R. R. Du, J. A. Simmons and J. L. Reno, 
Phys. Rev. Lett. 89, 076801 (2002); A. A. Bykov, J. Zhang, S. 



Vitkalov, A. K. Kalagin and A. K. Bakarov, Phys. Rev. B 72, 

245307 (2005); W. Zhang, H.-S. Chiang, M. A. Zudov, L. N. 

Pfeiffer and K. W. West, ibid 75, 041304(R) (2007). 
[4] A. T. Hatke, M. A. Zudov, L. N. Pfeiffer and K. W. West, Phys. 

Rev. B 79, 161308(R) (2009). 
[5] R. G. Mani at al, Nature 420, 646 (2002); M. A. Zudov, R. R. 

Du, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 90, 046807 

(2003) ; ibid 96, 236804 (2006); C. L. Yang at al, ibid 91, 
096803 (2003); A. V. Andreev, I. L. Aleiner, and A. J. Millis, 
ibid 91, 056803 (2003); A. Auerbach, I. Finkler, B. I. Halperin, 
and A. Yacoby, ibid 94, 196801 (2005); I. G. Finkler and B. I. 
Halperin, Phys. Rev. B 79, 085315 (2009). 

[6] W. Zhang, M. A. Zudov, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. Lett. 98, 106804 (2007); A. A. Bykov, J.-Q. Zhang, 
S. Vitkalov, A. K. Kalagin, and A. K. Bakarov, ibid 99, 1 16801 
(2007); W. Zhang, M. A. Zudov, L. N. Pfeiffer, and K. W. West, 
ibid 100, 036805 (2008). 

[7] V. I. Ryzhii, Sov. Phys. Solid State 11, 2078 (1970); V. I. 
Ryzhii, R. A. Suris, and B. S. Shchamkhalova, Sov. Phys. 20, 
1299 (1986); A. C. Durst, S. Sachdev, N. Read, and S. M. 
Girvin, Phys. Rev. Lett. 91, 086803 (2003); I. A. Dmitriev, 
A. D. Mirlin, and D. G. Polyakov, Phys. Rev. B 75, 245320 
(2007). 

[8] M. G. Vavilov and I. L. Aleiner, Phys. Rev. B 69, 035303 

(2004) . 

[9] M. G. Vavilov, I. L. Aleiner, and L. I. Glazman, Phys. Rev. B 
76, 115331 (2007). 
[10] I. A. Dmitriev, A. D. Mirlin, and D. G. Polyakov, Phys. Rev. 
Lett. 91, 226802 (2003); I. A. Dmitriev, M. G. Vavilov, I. L. 
Aleiner, A. D. Mirlin, and D. G. Polyakov, Phys. Rev. B 71, 
115316 (2005). 

[11] A. T. Hatke, H.-S. Chiang, M. A. Zudov, L. N. Pfeiffer, and 
K. W. West, Phys. Rev. B 77, 201304(R) (2008). 

[12] M. Khodas and M. G. Vavilov, Phys. Rev. B 78, 2453 19 (2008). 

[13] S. I. Dorozhkin, JETP Lett. 77, 577 (2003); S. I. Dorozhkin, 
J. H. Smet, K. von Klitzing, L. N. Pfeiffer, and K. W. West, 
ibid 86, 543 (2007); M. A. Zudov, Phys. Rev. B 69, 041304(R) 
(2004); M. A. Zudov, R. R. Du, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. B 73, 041303(R) (2006); R. L. Willett, L. N. Pfeif- 
fer, and K. W. West, Phys. Rev. Lett. 93, 026804 (2004); A. T. 
Hatke, H. -S. Chiang, M. A. Zudov, L. N. Pfeiffer, and K. W. 
West, ibid 101, 246811 (2008). 

[14] X. L. Lei and S. Y. Liu, Appl. Phys. Lett. 89, 182117 (2006); 
X. L. Lei, Phys. Rev. B 79, 115308 (2009). 

[15] I. A. Dmitriev, A. D. Mirlin, and D. G. Polyakov, Phys. Rev. 
Lett. 99, 206805 (2007). 

[16] I. V. Pechenezhskii, S. I. Dorozhkin, and I. A. Dmitriev, JETP 
Lett. 85, 86 (2007). 

[17] K. W. Chiu, T. K. Lee, and J. J. Quinn, Surf. Sci. 58, 182 (1976). 

[18] P. K. Tien and J. P. Gordon, Phys. Rev. 129, 647 (1963). 

[19] M. Moskalets and M. Buttiker, Phys. Rev. B 69, 205316 (2004). 



